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1.  Introduction 


During  the  past  decade  our  group  has  been  investigating  the  manner 
in  which  the  presence  of  water  in  the  atmosphere  affects  the  large-scale 
atmospheric  circulation.  Our  work  has  been  performed  under  a  number  of 
consecutive  contracts  with  the  Air  Force  Geophysics  Laboratory.  Our  more 
recent  work  has  centered  about  the  formulation  of  a  low-order  model  in 
which  water  occupies  a  key  postion.  We  have  constructed  the  model  so 
that  separate  properties  of  water  may  be  suppressed,  individually  or  in 
combination,  and  the  resulting  changes  in  the  circulation  may  be 
ascertained . 

The  model  is  an  extension  of  the  familiar  dry  two-layer 
quasi-geostrophic  beta-plane  model,  as  originally  formulated  by  Phillips 
(1951).  Isobaric  height  (or  stream  function)  and  temperature  at  any 
reference  level  may  be  chosen  as  the  basic  prognostic  dependent  variables 
in  the  dry  model.  Velocity  potential  and  individual  pressure  change 
appear  as  auxiliary  variables,  determined  diagnostically  from  the  height 
and  temperature  fields.  To  these  we  have  appended  total  dew  point  and 
ocean-surface  temperature  as  basic  prognostic  variables,  to  produce  a 
moist  model.  Ordinary  dew  point,  relative  humidity,  water-vapor  mixing 
ratio,  liquid  water  mixing  ratio,  and  saturation  mixing  ratios  at 
atmospheric  and  oceanic  temperatures  enter  as  additional  diagnostically 
determined  variables.  For  a  reference  level  we  have  chosen  the  earth's 
surface  (1000  mb). 
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The  diabatic  processes  include  viscous  drag,  sensible-heat  transfer, 
evaporation  and  precipitation,  and  short-wave  and  long-wave  radiation. 
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We  have  made  the  model  "low-order"  by  expressing  each  basic-variable 
field  as  a  double  Fourier  series,  and  then  truncating  each  series  to 
include  only  seven  terms,  just  as  we  previously  did  in  formulating  a 
low-order  dry  model  used  to  study  vacillation  (Lorenz,  1963).  A  complete 
description  of  the  model  appears  in  the  final  report  of  the  most  recently 
completed  contract  (F  1 9628-81 -K-001 ) ,  entitled  "A  low-order  model  of  a 
moist  general  circulation:  formulation  and  testing".  A  similar 
description,  accompanied  by  some  "synoptic  maps"  produced  by  the  model, 
appears  in  the  Journal  of  the  Atmospheric  Sciences  (Lorenz,  1984),  under 
the  title  "Formulation  of  a  low-order  model  of  a  moist  general 
circulation". 

In  that  article  we  describe  a  set  of  three  preliminary  numerical 
runs  made  with  the  model,  each  lasting  500  days.  In  the  initial  run  the 
chosen  cross-latitude  contrast  in  solar  heating  produces  a  zonal  flow, 
which  becomes  unstable  with  respect  to  Rossby-wave  perturbations  upon 
acquiring  sufficient  strength;  the  perturbations  subsequently  grow  to 
become  fully  developed  waves.  In  the  final  state  of  the  final  run,  which 
seems  to  be  representative,  the  surface  map  contains  a  succession  of 
cylones  and  anticyclones.  The  highest  relative  humidities  occurred 
somewhat  to  the  NNE  of  the  cyclone  centers,  while  the  maximum  rates  of 
precipitation  occur  somewhat  to  the  SE.  The  strongest  evaporation  occurs 
in  the  southern  portions  of  the  anticyclones. 

On  the  negative  side  of  the  ledger,  the  cyclones  and  anticyclones 
are  much  weaker  than  typical  systems  in  the  real  atmosphere.  On 
occasions  during  the  runs  the  model  systems  attain  more  realistic 
amplitudes.  We  concluded  that  the  model  behaved  in  a  qualitatively 
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reasonable  manner,  and  Chat,  with  possible  adjustments  in  the  numerical 
values  of  some  of  the  parameters,  it  was  ready  for  production  runs. 

Most  of  our  computations,  including  those  which  formed  the  three 
preliminary  runs,  have  been  carried  out  with  a  relatively  slow  system. 
This  report  describes  a  limited  number  of  runs  which  were  performed  on 
the  Cray-1  computer  at  the  National  Center  for  Atmospheric  Research. 

This  system  is  about  three  orders  of  magnitude  faster  than  the  system 
which  we  have  customarily  used,  and  it  will  integrate  the  equations  of 
our  model  for  one  year  of  simulated  time  in  about  1.3  seconds. 

In  the  sections  which  follow  we  first  describe  some  experiments 
using  various  combinations  of  values  of  the  model  parameters,  designed 
primarily  to  determine  suitable  combinations  for  use  in  subsequent 
experiments.  We  next  describe  a  very  long  run  performed  with  "suitable" 
parameter  values,  which  produced  fluctuations  with  strikingly  long 
periods.  Finally  we  present  some  comparison  runs  in  which  the  influence 
of  cloudiness  on  the  earth's  albedo  has  been  suppressed. 

2.  Some  adjustable  model  parameters 

As  we  have  already  noted,  the  cyclones  and  anticyclones  which 
develop  in  the  preliminary  runs  with  the  model  are  much  weaker  than  their 
real  atmospheric  counterparts.  In  our  detailed  write-up  we  mentioned 
that  this  discrepancy  might  have  resulted  from  an  excessively  high  static 
stability,  a  weak  cross-latitude  temperature  contrast,  or  some  other 
unidentified  cause.  The  experiments  to  be  described  in  this  section 
represent  a  first  attempt  to  account  for  the  discrepancy,  and  to 
eliminate  it  by  choosing  new  parameter  values,  i.e.,  to  tune  the  model. 
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The  model  has  a  preassigned  constant  lapse  rate  of  temperature, 
which  is  introduced  by  letting  the  temperature  in  each  vertical  column 
vary  as  p^f  where  p  is  pressure.  When  X  =  <  =  2/7,  the  lapse  rate  is 
dry  adiabatic.  In  the  preliminary  runs  we  set  X  =  0.175,  or  about  0.61 
times  the  dry-adiabatic  value.  We  noted,  incidentally,  that  the  use  of  a 
fixed  static  stability  was  one  of  the  main  shortcomings  of  the  model, 
which  we  hoped  eventually  to  remove,  at  least  to  the  extent  of  letting  X 
vary  with  time. 

The  model  atmosphere  extends  horizontally  over  an  infinite  channel, 
bounded  on  the  south  and  north  by  frictionless  walls  at  y  =  0  and  y  =  tt. 
The  basic  variables,  but  not  the  auxiliary  variables,  assume  the  form 
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X  -  Z  xn  «J>n  ,  (1) 
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Rossby  waves.  A  superposition  of  waves  on  a  suitable  zonal  flow  produces 
cyclones  and  anticyclones.  The  solar  heating  is  expressed  in  terms  of  a 
"planetary"  temperature  field  T*,  whose  value  at  any  location  is  the 
temperature  which  a  black  body  covering  the  entire  sky  would  have  to  be 
equivalent  to  the  sun.  In  the  preliminary  runs  we  chose 

T*  *  T o*  ®o  +  T3*  $3,  (3) 

so  that  T3*  measured  the  heating  contrast,  and  we  let  Tq*  =  273.0  and 
T3*  *  .0  (we  quote  all  values  of  T,  W,  and  S  in  degrees  Kelvin). 

A  feature  of  the  model's  behavior  which  became  apparent  in  the 
preliminary  runs,  and  must  be  recognized  if  the  tuning  process  is  to 
succeed,  is  the  high  sensitivity  of  the  temperature  T  to  the  planetary 
temperature  T*,  and  particularly  of  the  globally  averaged  surface 
temperature  T0  to  the  globally  averaged  planetary 

temperature  Tq*.  This  appears  to  result  from  a  cloud-albedo  feedback 
process.  Increased  solar  heating  not  only  heats  the  air,  but  also  tends 
to  lower  the  relative  humidity  and  the  accompanying  cloudiness  (although 
the  specific  humidity  increases),  thus  allowing  more  solar  energy  to 
reach  the  earth,  and  warming  the  atmosphere  still  more.  An  accompanying 
feature  is  the  very  slow  rate  at  which  Tq  approaches  its  ultimate  value 
from  a  prechosen  initial  state.  It  is  as  though  the  temperature  were 
attempting  to  reach  an  equilibrium,  and,  having  nearly  reached  it, 
discovered  that  the  rules  had  been  changed  and  that  it  must  approach 
another  equilibrium. 

One  property  of  the  real  atmosphere-ocean  system  which,  if  present 
in  the  model,  would  make  the  approach  to  equilibrium  very  slow  is  the 
high  ratio  R  of  the  heat  capacity  of  the  upper  layer  of  the  ocean  to  that 
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of  the  atmosphere.  For  an  oceanic  mixed  layer  70  meters  deep, 
interacting  with  the  atmosphere,  R  would  be  about  30.  In  our  preliminary 
runs,  in  order  to  reduce  the  required  computation  time,  we  let  R  =  1,  so 
that  the  ocean  might  more  properly  have  been  considered  to  be  wet  land. 
Despite  the  small  value  of  R,  the  approach  to  equilibrium  was  slow. 

The  runs  to  be  described  make  use  of  various  values  of  A,  Ty*,  T 3*, 
and  R.  All  other  adjustable  model  parameters  retain  the  values  used  in 
the  preliminary  runs. 

3.  Tuning  the  model 

Our  first  run,  which  we  shall  call  Run  1,  was  performed  to  determine 
whether  the  final  preliminary  run  had  indeed  reached  equilibrium,  and  it 
therefore  used  the  same  parameter  values.  Not  knowing  how  much  time 
might  really  be  needed  to  attain  equilibrium,  we  let  the  model  run  for  23 
years,  starting  with  a  small  perturbation  (  =  0,  T^  =  Wj  =  Sj  =  1.0) 

superposed  on  the  zonally  symmetric  (but  not  steady)  fields  ip  =  0, 

T  =  W  =  S  =  T*.  Fig.  1  shows  instantaneous  values  of  Ty  at  three-month 
intervals.  After  .1  brief  initial  drop,  Ty  increases  continually, 
although  somewhat  erratically,  attaining  a  peak  value  of  280.1  in  16.5 
years.  Thereafter  it  drops  slightly,  eventually  undergoing  fluctuations, 
witli  a  range  of  about  0.3,  about  a  mean  of  279.7. 

The  final  preliminary  run  and  the  new  run  are  entirely  consistent, 
except  that  Ty,  Wy,  and  Sy  are  about  one  degree  too  low  in  the  former. 

The  cross -la t i tude  temperature  contrast  T^  is  a^out  twice  T^*,  as  a 
result  of  the  sensitivity  of  T  to  T*,  while  the  Rossby  waves  are 
uniformly  weak. 
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We  next  performed  two  25-year  runs,  again  with  A  =  0.175  and 
R  =  1.0,  but  with  T3*  =  15.0,  and  with  T0*  =  273.0  in  Run  2  and  270.0  in 
Run  3.  Both  runs  finally  reached  quasi-equilibrium,  but  with  ’r(J 
fluctuating  around  284.9  in  the  former  run  and  271.1  in  the  latter.  Thus 
the  two  values  of  Tq  differed  by  more  than  4  times  the  corresponding 
difference  in  Tq*.  Moreover,  increasing  T 3*  from  5.0  to  15.0  increased 
To  by  more  than  five  degrees,  even  though  T0*  remained  fixed,  but  T3 
increased  only  slightly,  from  10.7  to  13.1.  Evidently  these  values  of  T 3 
are  comparable  to  the  minimum  value  needed  for  baroclinic  instability. 

The  waves  approximately  doubled  in  amplitude,  but  were  still  rather 
weak.  In  all  cases  the  waves  in  the  W-field  were  stronger  than  these  in 
the  T-field,  while  those  in  the  S-field  were  barely  detectable. 

The  new  runs  did  not  acquire  steady  states,  since  the  Rossby  w.cves 
continued  to  progress  eastward.  Moreover,  even  the  zonally  symmetric 
components  Tq f  T3 ,  and  T^  continued  to  fluctuate.  The  nature  of  these 
fluctuations  is  revealed  by  Fig.  2,  which  shows  ins tantaneous  values  of 
T(j  in  Run  2,  at  3-hour  intervals  for  25  days.  Accompanying  these  values 
are  horizontal  line  segments  indicating  the  times  when  Tj  and  T  4,  the 
cosine  phases  of  the  two  wave  inodes  in  the  T  field,  are  positive.  There 
are  periodic  oscillations  with  a  range  of  about  0.4  K  and  a  period  of 
about  11.2  days.  This  is  evidently  the  time  required  for  the  second  wave 
mode  ( Ti,  <61,  +  T5  $5),  traveling  at  about  20  m  s_  to  move  one  wave 
length  (5750  km)  farther  than  the  first  wave  mode(Tj  4>j  ♦  T^  $2)* 
traveling  at  about  14  m  s-^.  These  modes  interact  with  T3  through  an 
advective  (quadratic)  non  1  i near i  ly ,  while  Tq  affects  Tq  through  a 
higher-degree  nonlinearity.  Thus  vacillation  is  occurring. 
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Superposed  on  this  oscillation  is  a  weaker  fluctuation  with  a  period 
of  about  1.1  days,  which  appears  to  be  the  time  required  for  the  second 
wave  mode  to  travel  1/3  wave  length.  This  time  is  significant  because, 
in  integrating  the  equations,  we  transform  back  and  forth  between  Fourier 
modes  and  grid  points  during  each  time  step.  The  west-east  spacing  of 
the  grid  points  is  three  per  wave  length.  We  evaluate  the  high-degree 
nonlinear  terms  which  enter  the  radiation  and  moist  thermodynamic 
formulas  at  grid  points;  the  total  effect  of  these  terms  on  T q  differs 
according  to  whether  a  wave  trough  (or  ridge)  is  located  close  to  a  grid 
point  or  almost  midway  between  grid  points. 

We  could  virtually  eliminate  this  effect  by  using  a  denser  network 
of  grid  points,  but  this  would  conflict  with  our  aim  of  keeping  the  model 
as  small  as  possible.  We  prefer  to  keep  the  model  as  it  is,  and  regard 
it  as  possessing  weak  periodic  topography  with  a  period  of  one  grid 
interval.  The  effect  of  such  topography  would  be  to  produce  oscillations 
like  the  shorter-per i od  ones  in  Fig.  2. 

Similar  doubly  periodic  oscillations  occur  in  Runs  1  and  3.  The 
apparent  15-month  periodicity  occurring  during  the  last  6  years  in 
Fig.  1  is  actually  a  short-period  oscillation,  which  has  been  aliased  by 
being  sampled  only  once  every  three  months.  The  two  short  periods  in 
Fig.  2  appear  incommensurable,  but  there  is  no  evidence  of  aperi odi c i ty . 

For  Run  4  we  have  increased  A  to  0.2,  while  letting  T y  =  273.0, 

Tj  =  10.0,  and  R  =  1.0.  In  Fig.  3  the  vertical  line  segments  show  the 
total  range  of  T(j  during  each  three-month  inLerval.  Again  T  ^  appears  to 
be  seeking  an  equilibrium,  but  the  fluctuations  are  more  pronounced  than 
those  appearing  in  Fig.  1.  Moreover,  they  are  not  merely  aliased 
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On  the  basis  of  Runs  1-6  we  have  decided  that  the  values  A  =  0.2, 

Ty*  =  273.0,  T 3*  =  10.0,  and  R  =  10.0  are  suitable  for  further  runs.  It 
appears  unlikely  that  such  runs  will  ever  settle  down  to  steady  behavior, 
so  that  a  considerable  variety  of  initial  states  should  be  compatible 
with  the  climate.  These  should  include  any  state  in  the  later  years  of 
Run  6,  but  suitable  states  also  appear  to  be  reached  after  a  few  months 
from  the  simple  initial  state  consisting  of  a  small  perturbation, 

=  1.0,  superposed  on  the  zonally  symmetric  state  where 
vo  =  0,  To  =  287.0,  W0  =  284.5  ,  S0  =  286.5,  1(13  =  0,  T3  =  6.0,  W3  =  4.5, 

=  5.5,  and  =  Tb  =  Wb  =  Se  =  0. 

4.  Low-frequency  fluctuations 

The  later  years  of  Run  4  reveal  considerable  irregular  variability, 
with  peak  values  of  Ty  occurring  at  1-  to  2-year  intervals.  In  Run  6, 
where  R  has  been  increased  to  10.0,  peaks  show  some  sign  of  occurring  at 
5-  to  10-year  intervals,  and  a  general  downward  ^ -  nd  is  evident.  To 
determine  whether  the  low-frequency  fluctuations  and  the  downward  trend 
continue,  we  have  extended  Run  6  for  370  additional  years;  the  entire  450 
years  constitute  Run  7. 

Fig.  6  shows  tlie  successive  annual  averages  of  Ty  for  the  final  400 
years  of  Run  7,  which  begin  shortly  after  the  peak  value  in  Fig.  5  has 
been  atttained.  The  downward  trend  soon  ceases,  and  there  appear  to  be 
irregular  variations  on  all  time  scales.  The  total  range  of  Ty  exceeds 
two  degrees.  Changes  from  one  year  to  the  next  seldom  exceed  0.4,  while 
increases  or  decreases  of  l .0  or  more  generally  require  10  years  or 
longer.  Ten-year  averages  vary  nearly  as  much  as  annual  averages,  while 
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short-period  oscillations.  Such  osci llations,  with  a  range  of  about 
0.5  K,  could  possibly  account  for  the  lengths  of  the  segments,  but  the 
presence  of  completely  non-overlapping  segments,  some  less  than  a  year 
apart,  shows  that  longer-period  fluctuations  are  also  present.  Similarly 
drawn  segments  for  the  final  years  of  Run  1  would  show  no  variability  at 
ail. 

Meanwhile  T3  fluctuates  about  a  mean  value  of  about  5.6;  this  value, 
about  half  of  that  in  Runs  1-3,  is  apparently  large  enough  to  produce 
baroclinic  instability  with  the  steeper  lapse  rate  (A  =  0.2).  The  waves 
are  now  about  4  times  as  strong  as  in  Run  1,  and  have  entered  the 
atmospheric  range. 

To  observe  the  effect  of  using  a  more  realistic  oceanic  heat 
capacity,  we  have  performed  Run  5,  which  is  like  Run  4  except  that 
R  =  25.0.  The  result  is  indicated  by  the  dots  in  Fig.  4,  which,  as  in 
Fig.  1,  show  instantaneous  values  of  Tq.  After  15  years,  Tq  seems  no 
closer  to  equilibrium  than  it  was  after  2  years  in  Run  4.  Determining 
the  final  equilibrium  could  be  a  laborious  task. 

As  a  compromise  we  have  performed  Run  6,  with  R  =  10.0  and  with  the 
other  parameters  unchanged  from  Run  4,  and  have  used  the  final  conditions 
of  Run  4  as  initial  conditions.  In  Fig.  5  the  vertical  segments  indicate 
the  total  range  of  Tq  in  each  of  80  successive  years.  Within  each  year 
the  annual  range  is  comparable  to  1.0.  The  annual  mean  increases,  with 
one  lapse,  for  almost  20  years;  thereafter  it  remains  reasonably  steady 
for  another  20  years.  During  the  later  years  there  is  a  general  downward 
trend,  and  Tq  loses  a  full  degree.  The  annual  mean  values  of  Wq  behave 
similarly,  and  in  fact  almost  exactly  fit  the  linear  relation 
T0-287 . 2  -  1.2(Wq-284.5) . 
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even  50-year  averages  are  insufficient  for  estimating  true  climatic 
means;  during  years  100-150  the  average  Tq  is  287.1,  while  during  years 
200-250  it  is  287.8.  Apparently  the  variations  of  Tq  are  dominated  by 
very  low  frequencies. 

A  one-degree  change  in  a  local  temperature  is  not  a  very  striking 
event.  The  significance  of  the  changes  shown  in  Fig.  6  is  that  they  are 
changes  in  the  globally  and  annually  averaged  temperature,  and  that 
variations  of  the  global  annual  average  temperature  of  the  real 
atmosphere  seem  to  be  of  comparable  magnitude,  occurring  on  comparable 
time  scales  (cf.  Pan  and  Oort,  1983).  In  all  likelihood  a  warming  or 
cooling  of  the  real  atmosphere,  proceeding  more  or  less  continually  for 
10  or  20  years,  and  totaling  one  degree,  would  be  interpreted  by  many 
observers  as  a  climatic  change,  and  attempts  might  be  made  to  determine  a 
cause. 

In  the  model  there  is  no  special  cause  for  these  changes;  they 
simply  represent  the  model's  natural  variability.  There  are  no  changes 
whatever  in  the  solar  output,  nor  in  any  other  conditions  external  to  the 
atmosphere  and  ocean.  While  a  model  as  simple  as  this  one  can  yield  only 
limited  information  about  the  real  atmosphere,  it  at  least  shows  that 
changes  having  the  magnitude  and  time  scale  of  those  encountered  in  Run  7 
do  not  necessarily  demand  external  causes. 

5.  Comparison  runs 

One  of  our  principal  motives  tor  constructing  the  low-order  moist 
model  was  our  wish  to  investigate  the  influence  of  water  on  the 
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circulation  of  the  atmosphere,  by  comparing  ordinary  runs  with  special 
runs  in  which  the  influence  of  water  is  treated  differently.  Here  we 
shall  compare  the  runs  described  in  the  previous  section  with  runs  in 
which  the  influence  of  clouds  on  the  earth's  albedo  has  been  suppressed. 
In  the  new  runs  we  have  given  the  albedo  the  constant  value  of  0.297; 
this  is  the  value  which  it  will  assume  in  an  ordinary  run  when  T  =  285.0. 

In  the  first  three  new  runs  we  have  let  Tq*  =  270.0,  273.0,  and 
276.0,  while,  as  in  Run  4,  A  =  0.2,  T  3*  =  10.0,  and  R  =  1.0.  In  each  run 
quas i -equi li br i um  is  attained  in  less  than  two  years.  The  equilibrium 
values  of  Tq  are  about  281.0,  284.1,  and  287.9,  while  the  values  of  T3*, 
and  the  intensities  of  the  cyclones  and  anticyclones,  are  comparable  to 
those  in  Run  4.  The  high  sensitivity  of  T0  to  T  q*  is  no  longer  present. 
Clearly  this  sensitivity  and  the  slow  approach  to  equilibrium  in  the 
ordinary  runs  are  associated  in  some  way  with  the  variable  albedo. 

In  the  next  three  new  runs  we  have  let  T3*  =  4.0,  8.0,  and  12.0, 
while,  again  as  in  Run  4,  A  =  0.2,  Tq*  *  273.0,  and  R  =  1.0.  The  first 
new  run  settles  down  after  slightly  more  than  a  year  to  steady  zonally 
symmetric  equilibrium,  with  Tq  =  283.0  and  T3  =  4.9;  evidently  this  value 
of  T3  is  insufficient  for  baroclinic  instability.  In  the  other  new  runs 
tlie  waves  persist,  becoming  somewhat  weaker  than  those  in  Run  4  in  the 
run  where  T3*  =  8.0,  and  acquiring  comparable  intensity  when  T3*  =  12.0. 
The  respective  values  about  which  T3  continues  to  fluctuate  are  about  5.2 
and  5.5,  which  evidently  exceed  the  baroc 1 i ni c-ins tabi li ty  criterion.  It 
thus  appears  that,  in  the  present  model,  a  cloud-albedo  process  can 
enhance  the  ability  of  a  weak  cross-latitude  heating  contrast  to  produce 
baroclinic  instability. 
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In  our  final  new  run  we  have  increased  R  to  10.0,  while  again 
letting  A  =  0.2,  Tq  =  273.0,  and  T 3*  =  10.0.  With  this  high  oceanic  heat 
capacity  nearly  10  years  are  needed  for  Tq  to  attain  quasi -equi librium, 
but  the  approach  is  still  considerably  more  rapid  than  in  Run  4,  where 
the  heat  capacity  is  low  (R  =  1.0)  but  the  cloud-albedo  feedback  is 
present . 

Additional  alterations  of  physical  processes  which  might  be  made  in 
comparison  runs  are  almost  limitless.  The  examples  which  we  have 
described  are  sufficient  to  show  that  simple  changes  can  have  pronounced 
effects. 

6.  Concluding  remarks 

Provided  that  one  possesses  a  numerical  model  to  work  with,  one  can 
test  the  importance  of  a  physical  process  which  is  present  in  the  model 
by  determining  the  typical  behavior  of  the  model,  and  subsequently 
inactivating  the  physical  process  in  question  and  observing  how  the 
behavior  of  the  model  has  been  altered.  We  have  demonstrated  that  this 
technique  can  be  used  even  with  a  very-low-order  model.  Our  model  was 
designed  to  deal  with  an  atmosphere  which  possesses  an  important 
constituent  which  occurs  with  variable  concentration  and  in  more  than  one 
phase.  In  our  applications  we  have  identified  the  atmosphere  with  the 
earth's  atmosphere,  and  the  constituent  with  water. 

We  have  presented  one  example  indicating  that  the  equilibrium  state 
about  which  the  model  atmosphere  will  fluctuate,  and  the  rate  at  which  it 
will  approach  this  equilibrium  from  an  unbalanced  state,  are  strongly 


influenced  by  a  cloud-albedo  feedback  process.  This  result  may  not  hold 
up  too  well  in  the  real  world,  since  it  appears  that  the  model  may 
exaggerate  either  the  influence  of  clouds  on  the  albedo  or  the  extent  to 
which  cloudiness  varies.  Nevertheless,  we  have  shown  that  the  model  may 
be  tuned  to  fit  the  earth's  atmosphere  more  closely.  Additional  tuning, 
perhaps  by  changing  the  value  of  the  parameter  representing  the  "half 
life"  of  liquid  water,  i.e.,  the  time  during  which  liquid  water  tends  to 
persist  as  clouds  before  falling  out  as  rain,  could  bring  about  a  still 
better  fit.  Also  there  are  relevant  physical  processes  which  for 
simplicity  we  have  omitted,  but  which  could  easily  be  introduced. 

Beyond  these  considerations,  we  have  devoted  considerable  attention 
to  one  of  the  more  interesting  and  unexpected  model  phenomena  to  be 
discovered,  namely,  the  pronounced  long-term  fluctuations  of  the  globally 
averaged  temperature,  which  resemble  climatic  changes.  A  natural 
conclusion  is  that  the  model  might  find  use  as  a  climate  model.  It  would 
not  be  too  difficult  to  change  to  spherical  geometry,  and  introduce  such 
features  as  seasonal  variations  in  the  solar  forcing,  and  the  presence  of 
polar  ice,  which  a  climate  model  ought  to  possess.  The  great  advantage 
of  our  model  is  its  speed;  with  a  large  computer,  runs  of  1,000  or  even 
10,000  years  are  entirely  feasible.  We  venture  a  guess  that  the  model, 
with  suitable  additions,  could  prove  capable  of  creating  an  ice  age  from 
an  interglacial  initial  state,  or  vice  versa,  even  though  the  numerical 
integration  must  proceed  in  1.5-hour  time  steps. 


15 


i 


; . 

.**-.*Vk**_»" 


References 

Lorenz,  E.  N.,  1963:  The  mechanics  of  vacillation.  -J.  Atmos.  Sci., 
20,  448-464. 

Lorenz,  E.  N.,  1984:  Formulation  of  a  low-order  model  of  a  moist 
general  circulation.  J,  Atmos.  Sci.,  41 ,  1933-1945. 

Pan,  Y.  H.  and  A.  H.  Oort,  1983:  Global  climate  variations 
connected  with  sea  surface  temperature  anomalies  in  the  eastern 
equatorial  Pacific  Ocean  for  the  1958-73  period.  Mon.  Wea.  Rev.,  Ill, 
1244-1258. 

Phillips,  N.  A.,  1951:  A  simple  three-dimensional  model  for  the 
study  of  large-scale  extratropical  flow  patterns.  J.  Meteor.,  f$, 

381-394. 


•• 


0 


Fig.  4.  Instantaneous  values  of  Tq  at  3-month  intervals  in  Run  5,  where 
X  =  0.2,  Tq*  =  273.0,  T3*  =  10.0,  and  R  =  25.0. 
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